Generation of nonclassical photon states using a superconducting qubit in a microcavity 
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Based on the interaction between the radiation field and a superconductor, we propose a way to engineer 
quantum states using a SQUID charge qubit inside a microcavity. This device can act as a deterministic single 
photon source as well as generate any Fock states and an arbitrary superposition of Fock states for the cavity 
field. The controllable interaction between the cavity field and the qubit can be realized by the tunable gate 
voltage and classical magnetic field applied to the SQUID. 

PACS numbers: 42.50.Dv, 74.50.+r, 42.50.Ct 



The generation of quantum states of the radiation field 
has been a topic of growing interest in recent years. This 
is because of possible applications in quantum communica- 
tion and information processing, such as quantum networks, 
secure quantum communications, and quantum cryptogra- 
phy [l|]. Based on the interaction between the radiation field 
and atoms, many theoretical schemes have been proposed for 
the generation of Fock states and their arbitrary super- 

positions 1 4, 5]. Experiments have generated single-photon 
states in quantum dots @], atoms inside a microcavity Q], 
and other systems fH]. A superposition of the vacuum and 
one-photon states can also be experimentally created by trun- 
cating an input coherent state or using cavity quantum electro- 
dynamics 1 9]. However, how to generate an arbitrary photon 
state by virtue of the interaction between the radiation field 
and solid state quantum devices seems to be unknown both 
theoretically and experimentally. Recent progress in super- 
conducting quantum devices (e.g., 11(1 ll 110 makes it possi- 
ble to do quantum state engineering experiments in these sys- 
tems, and also there have been proposals on superconduct- 
ing qubits interacti ng w ith the nonclassical electromagnetic 

field EESHHtiEO 

Here, we present an experimentally feasible scheme to gen- 
erate quantum states of a single-mode cavity field in the mi- 
crowave regime by using the photon transition between the 
ground and first excited states of a macroscopic two-level sys- 
tem formed by a superconducting quantum interference de- 
vice (SQUID). This artificial two-level "atom" can be easily 
controlled by an applied gate voltage V g and the flux <i> c gener- 
ated by the classical magnetic field through the SQUID (e.g., 
III4LII8II ). The process of generating photon states in this de- 
vice includes three main steps: (i) The artificial atom operates 
at the degeneracy point by choosing appropriate values for V g 
and $ c . There is no interaction between the quantized cavity 
field and "atom" at this stage, (ii) Afterwards new V g and $ c 
are selected such that the cavity field interacts resonantly with 
the "atom" and evolves during a designated time, (iii) The 
above two steps can be repeated until a desired state is ob- 
tained. Finally, the flux $ c can be adjusted to a special value, 
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FIG. 1: (a) A charge qubit formed by a SQUID device, equivalent 
to a controllable macroscopic two-level system, is placed into a su- 
perconducting microwave cavity in (b). The coupling between the 
quantized cavity field and qubit system is realized via the magnetic 
flux $x through the SQUID. 



then the interaction is switched off, and the desired photon 
state appears in the cavity. This process is similar to that of a 
micromaser [2] and it is described below. 

Model. — The macroscopic two-level system studied here is 
shown in Fig.^Ca). A SQUID-type superconducting box with 
n c excess Cooper-pair charges is connected to a superconduct- 
ing loop via two identical Josephson junctions with capacitors 
Cj and coupling energies Ej . A controllable gate voltage V g 
is coupled to the box via a gate capacitor C g . We assume 
that the superconducting energy gap A is the largest energy. 
Then, at low temperatures, the quasi-particle tunneling is sup- 
pressed and no quasi-particle excitation can be found on the 
island. Only Cooper pairs coherently tunnel in the supercon- 
ducting junctions. The above assumptions are consistent with 
most experiments on charge qubits. Then the standard Hamil- 
tonian 1 18Q is 



H qh = 4£ ch (n c 



2£jcos( — ^ ) cos9,(l) 



where <£>x is the total flux through the SQUID loop and <P 
the flux quantum. Thus, the superconducting loop is used 
to control the Josephson coupling energy by adjusting the 
flux through this loop. Below, we show that it can also 
switch on and off the qubit-field interaction. The dimension- 
less gate charge, n g = C g V g /2e, is controlled by V g . The 



2 



single-electron charging energy is = e 2 /2(C g + 2Cj). 
= (0i + 02)/2 is the quantum mechanical conjugate of the 
number operator n c of the Cooper pairs on the box, where <^ 
(i = 1, 2) is the phase difference for each junction. The su- 
perconducting box is assumed to work in the charging regime 
with condition k B T < Ej < E cXl < A where T and k B 
are temperature and Boltzmann constant respectively. If the 
gate voltages is near a degeneracy point n g = 1/2, the super- 
conducting box is a charge qubit II 811 . which is a controllable 
two-level system characterized by the two lowest charge states 
\g) (for n c = 0) and \e) (for n c = 1). However, if the quasi- 
particle excitation cannot be completely suppressed, a contin- 
uum of low-lying quasi-particle states will be present, and the 
Hamiltonian cannot be reduced to a system with two en- 
ergy levels even when the gate voltage is near the degeneracy 
point ||19ll . 

Now we further consider that the qubit is placed in a 
single-mode microwave superconducting cavity, depicted in 
Fig-Efb), the flux $x through the SQUID can be expressed 
as ll3 E 01 $x = $ c + * q where the flux $ c and 
= r\a + rj* a' are generated by a classical applied mag- 
netic field and the quantized cavity field, respectively. Here 
rj = j s u(r) • ds and u(r) is the mode function of the cavity 
field, with annihilation (creation) operators a(aJ), and S is 
the surface defined by the contour of the SQUID. Considering 
the above, we obtain 



H = hwa^a + E z a z 
- E] (er + + cr_ ) cos 



(2) 



— ($ G 7 + r\a + rj* a)) 



where the first two terms represent the free Hamiltonians of 
the cavity field with frequency uo — AE^/h and the qubit 
with the energy E z = — 2E C ^(1 — 2n g ), / is the identity oper- 
ator. The third term is the nonlinear photon-qubit interaction 
which is switchable by the flux <1> C . The charge excited state 
|e) and ground state \g) correspond to the eigenstates | J.) and 
1 1) of the spin operator a z , respectively. The cosine in Eq. (0 
can be further decomposed into classical and quantized parts, 
and the quantized parts sin[7r(?7 a+H.c.)/ $o] an d cos[7r(?7 a+ 
H.c.) I $o] can be further expanded as a power series in a (cr). 
Here, the single photon transition between the states \e,n) 
and \g, n + 1) satisfies the condition (^\ri\/^o)V n + 1 <C 1> 
where n is the number of photons; therefore all higher or- 
ders of 7r|r7|/$o can be neglected and only a single-photon 
transition is kept in the expansion of Eq. (0. Using the no- 
tation for trapped ion systems (e.g., 12011 ), the first red (blue) 
sideband excitations (3a<7 + + H.c. (/3acr_ + H.c.) for inter- 
actions of the cavity field and the qubit fl3ll . with photon- 
qubit coupling constant /3 = (wriEj/^o) sin(7r$ c /$o)> can 
be obtained by adjusting the gate voltages V g and the flux $ c . 
They correspond to 2E Z = fiui (2E Z = ~hio) and dimension- 
less gate charge n g — 1 (n g = 0). Also £{0+ + a ) with 
£ = Ej cos(7r<i> c /$o) is called the carrier Ml 311 - which corre- 
sponds to n g = 1/2. The Hamiltonian (|2), with the above 
assumptions, is our model. 



Preparation process. — We choose |0, g) as our initial 
state, where the cavity field is in the vacuum state |0) and 
the qubit is in the ground state \g). The goal is to prepare an 
arbitrary pure state of the cavity field 
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(3) 



n=0 



n=0 



where |n) denotes the Fock states of the cavity field with ex- 
citation number n = 0, 1, 2, • • • . A Fock state \m) with m 
photons is a special case of Eq. Q with conditions c n = for 
all n ^ m with < m < N. 

Thermal photons in the cavity have to be suppressed in or- 
der to obtain the vacuum state |0). In the microwave region 
0.1 ~ 15 cm, the mean number of thermal photons (nth) sat- 
isfies 3.0 x 10~ 208 < (n t h) < 0.043 at T = 30 mK, and 
1.7 x lO" 104 < (n th ) < 0.26 at T = 60 mK. These tempera- 
tures can be obtained experimentally (e.g., in (DlEl)- 

After the system is initialized, two different processes are 
required to engineer the state of the cavity field. The first pro- 
cess involves rotating the qubit state, but keeping the cavity 
field state unchanged. This stage can be experimentally re- 
alized by tuning the gate voltage and classical magnetic field 
such that ?i g = 1/2 and $ c = 0; then the time evolution oper- 
ator Uc(t) of the qubit in the interaction picture is 

Uo(t) = cos(Oit)/ + iain(ni*)(|0)(e| + \e)(g\) (4) 

where fii = Ej/h. The subscript "C" in Uc(t) denotes the 
carrier process, which can superpose two levels of the qubit, 
and it can also flip the ground state \g) or excited state |e) to 
each other, after a time t = ir(2p — l)/2fii, with positive 
integer p. 

The second process is the first red (blue) sideband excita- 
tion, which can be realized by tuning the gate voltage and 
classical magnetic field such that n g = 1 (n g = 0) and 
$ c = <1>q/2. Thus, in the interaction picture, the time evo- 
lution operators U-R,{t) for the red (JJ&(t) for the blue) of the 
cavity field and qubit can be expressed 1 22] as 



U R (t) = R ee (t)\e)(e\+R gg (t)\9)(g\ 
- iR ge (t)\g){e\-iR eg (t)\e)(g\ 



Us(t) = Rgg 

- iR, 



(t)\e)(e\ + R ee (t)\g)(g\ 
ge (t)\e)(g\-iR eg (t)\g)(e\ 



(5) 



(6) 



with R eg (t) = e l9 sin Mfi 2 |waa+J /vaa) a, R ge (t) = 
Bt g (t), Ree{t) = cos{\n 2 \tVm^, and R gg {t) = 
cos ^| £^2 'at a^j , where we have assumed that fl2 = 

TrrjEj/h&o = I £^2 1 e* 6 ", in which the phase 9 depends on the 
mode function of the cavity field u(r). The red sideband exci- 
tation described by operator Unit) can entangle \g,n + 1) 
with |e,n), or flip \g,n + 1) to \e,n) and vice versa, by 
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choosing the duration of the interaction between the cavity 
field and the qubit. From Eq. (0, it is easy to verify that 
the emission probability P g of the upper level for the qubit 
is P g = sm 2 (\Q 2 \tVnTT). We find that P g = 1 when 
|^2^V 71 + 1 = 7r(2fe — l)/2, with positive integer k. So 
when t = 7r(2fc - l)/(2\Q 2 \\/n + 1), there are n + 1 pho- 
tons in the cavity and the qubit is in its ground state. The 
first blue sideband excitation, denoted by Ub it), can entangle 
state |e, n+1) with state \g, n), or flip |e, n + 1) to |g, n) and 
vice versa. Below we use the carrier and the first red sideband 
excitation, represented by {7c (i) and C/r (^), as an example 
showing the generation of an arbitrary quantum state of the 
cavity field. 

Using the quantum operations Uc(t) and Unit) in Eqs. (0} 
and (0, the single photon state |1) can be generated from 
the initial vacuum state |0). That is, we can first flip the 
ground state of the qubit to the excited state when the con- 
dition Qit\ = ir/2 is satisfied for the carrier Uc{t\), then 
we turn on the first red sideband excitation U R (t 2 ) and let the 
photon-qubit system evolve a time t 2 satisfying the condition 
| O2 ! *2 = t/2. Finally, we adjust the classical magnetic field 
such that $ c = 0; thus the interaction between the cavity field 
and qubit vanishes, and a single-photon state exists in the cav- 
ity, that is, 

\l)®\9)=U K (h\)Uc(h)\Q)®\g). (7) 

Also any Fock state \m) can be easily created from the vac- 
uum state |0) by alternatively turning on and off the quantum 
operations in Eqs. ( 1415 \ to excite the qubit and emit photons 
during the time interval T. The latter is divided by 2m subin- 
tervals n, r 2) • • • , T21-1, t% • ■ ■ , r 2m -i , T 2m which satisfy 
conditions |£2i|t2J-i = tt/2 and \£1 2 \t 2 i^/1 + 1 = tt/2 where 
I = 1, • • • , m. This process can be described as 

|m}<8> Iff) = U R (T 2m )Uc(T 2m -i) ■ ■ ■ U n (r 2 )Uc(n)\0) <g> \g). 

(8) 

Finally, the classical magnetic field is changed such that 3> c = 
$o> and an n-photon state is provided in the cavity. 

Our next goal is to prepare superpositions of different Fock 
states (e.g., ai|0) + 0:2 11)) for the vacuum |0) and single 
photon |1) states. This very important state can be deter- 
ministically generated by two steps, Uc(t'i) and U R (t' 2 ), with 
t' 2 = 7r/2|Q 2 |; that is 

(a x |0) + a a |l» <8> Iff) = U R (&)U c (£)\0) ® Iff) (9) 

where the operation time t[ determines the weights of the 
coefficients of the superposition a,\ — cos(Oit^) and a 2 = 
e sm(flit[). If the condition t[ = 7r/4f2i is satisfied, then 
we have a superposition (|0) + e~ l9 \l))/y/2 with equal prob- 
abilities for each component and the relative phase between 
them can be further specified by the phase of the mode of the 
cavity field. 

An arbitrary target state (0 can be generated from the ini- 
tial state by alternatively switching on and off the carrier and 
first red sideband excitation during the time X", which can be 



divided into 2n subintervals r{, • • ■ , r 2n . That is, the target 
state can be deterministically generated as follows 

JV 

\i>)=Y,Cn\n,9)=U(T')\0,g), (10) 

n=0 

where U(T') is determined by a sequence of time evolu- 
tion operators associated with chosen time subintervals as 

U(T>) = U R (T> n )Uc(T> n _ 1 )---U R (T>)U c (T[). Therefore, 
the coefficients c n are 

c n = (9,n\U R (T^)Uc(^ n _ 1 )---U R (r^Uc(ri)\0,g). 

(ID 

Reference J4J] has explicitly discussed how to adjust the 
rescaled times to obtain the expected state by solving the in- 
verse evolution of Eq. dlOi . Ideally, any state of the cavity 
field can be created according to our proposal by adjusting 
the gate voltage, classical magnetic field, and duration of the 
photon-qubit interaction. It is very easy to check that the state 
(0 can also be created by the carrier and blue sideband exci- 
tation whose time evolutions Uc(t) and U R (t) are described 
by Eqs. @ and (|6j. 

Environmental effects. — We now discuss the environmen- 
tal effects on the prepared states, which are actually limited by 
the following time scales: the relaxation time T\, the prepara- 
tion time r e of the excited state, and the dephasing time T 2 of 
the qubit, the lifetime r p of the photon and an effective interac- 
tion time ri™^ which corresponds to the transition from |n, e) 
and 1 71 + 1, ff). If Ti, t p r e , ri n \ then the Fock states 

(n) 

can be prepared. If the condition Tj., T 2 , t v ^> r e , Tc is 
satisfied, then the superposition can also be obtained. 

Now let us estimate the photon number of the obtainable 
Fock state in a full-wave cavity. In microwave experiments, 
it is possible to obtain very high-Q superconducting cavi- 
ties, with Q values around 3 x 10 8 to 5 x 10 10 jl Eill . 
which correspond to the lifetimes of the microwave region 
from 0.001 < t p < 0.15 s to 0.167 < r p < 25 s. The 
parameters of the charge qubit 12411 without the cavity are 
2Ej/h = 13.0 GHz (so the operation time corresponding to 
a completely excited qubit is about r e sw 3.8 x 10 -11 s). The 
lifetime of the excited-state for the qubit T\ = 1.3 x 10~ 6 s, 
i.e. T e <C T\. For an estimate of the interaction coupling 
between the cavity field and the qubit, we assume that the 
cavity mode function is taken as a standing-wave form such 
as B x — —iyjhjj 1 'eqV 'c 2 (a — a') cos(fcz), with polarization 
along the normal direction of the surface area of the SQUID, 
located at an antinode of the standing-wave mode; then the 
interaction between the cavity field and the qubit reaches its 
maximum and the interaction strength can be expressed as 
= n\r)\Ej/$> = {-KSEj/e&^ytfrw/euV. For example, 
if the wavelength of the cavity mode is taken as Ai = 0.1 
cm, then 7r|?7|/$ ~ 7.38 x 10~ 5 -C 1, where the dimen- 
sion of the SQUID is taken as 10 /im and the mode function 
u(r) is assumed to be independent of the integral area be- 
cause the dimension of the SQUID, 10 /im, is much less than 
0.1 cm, for the wavelength of the cavity mode. In this case, 
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FIG. 2: Ratio t^/t^, versus photon number n > 1, of the life- 
time Tp™' of the photon number state \n) and the effective operation 
time Tc ■ The latter corresponds to the transition from |n, e) to 
| ri. + 1, g) for a 10/im x 10/im SQUID in the full-wave cavity. 

Tc " 1 w 5.0 x 10~ 7 s, which is less than one order of magni- 
tude of the excited lifetime T\. This means that the qubit in its 
excited state can emit a photon before it relaxes to its ground 
state. But if we take the dimension of the SQUID as 1 /im, the 
coupling between the cavity field and the qubit is two orders 
of magnitude smaller than for the 10 /im SQUID, and then 
the interaction time is 5.0 x 10~ 5 > T\. Therefore, in this 
case, the qubit relaxes to the ground state before the photon 
can be emitted from the qubit, and thus it is difficult to obtain 
a photon state. In Fig. (0, we plot the ratio Tp^/Vc be- 
tween the lifetime lE5ll Tp 1 ^ — r p /n (here, Tp^ ~ r p ) of the 
Fock state |n) for the zero-temperature environment and the 
effective operation time of transitions from state \n, e) to 
\n+ 1, g) for different values of Q and for A = 0.1 cm. Fig. [2] 
shows that the photon number of the prepared Fock states can 
reach 10 2 in the above mentioned high-Q cavity. But if the 
Q values are less than 10 7 , it might be difficult to prepare a 
photon state with our estimated coupling. We also find that 
a longer microwave in the full-wave cavity corresponds to a 

(n) 

longer Tc ' for a fixed Q, which means that it is easy to cre- 
ate photon states for shorter microwaves. For example, if the 
wavelength is taken as A = 1 cm, then the coupling between 
the qubit-photon in the full-wave cavity might not be strong 
enough for generating Fock states within the currently known 
experimental data for T\. So for longer microwaves, we can 
make a smaller cavity and place the qubit where the qubit- 
photon interaction is maximum. 

If we want to prepare a superposition of different Fock 
states of the cavity field, we need to consider T2, which is 
of the order of a few ns (e.g., 5 ns in HHn . Then the survival 
time of the entangled state between the cavity field and qubit, 
which is required for the preparation of the superposition of 
the Fock states, maybe be very short. With the improvement 
of read-out techniques, a longer dephasing time can make our 
proposal far more realizable. 

Discussions. — We propose a scheme for deterministic ally 
generating nonclassical photon states via the interaction of 
photons and a charge qubit. Indeed, the Fock state can be 



prepared with current technology. The superposition would 
be easier to obtain by increasing the dephasing time T2 and 
the qubit-photon coupling strength. Our discussions above are 
based on experimental values for T\ and T 2 without the cav- 
ity; the decoherence may become shorter when the SQUID is 
placed inside the cavity. Further, in order to obtain a stronger 
coupling, the following steps would help to increase the qubit- 
field coupling strength: i) decrease the volume V of the cavity; 
ii) increase the area S of the SQUID; iii) increase the Joseph- 
son coupling energy Ej under the condition -C E c ^. 
We can also put a high permeability /1 material inside the 
SQUID loop I14II . then the qubit-field coupling strength can 
increase to because the relative permeability in ferro- 
magnetic materials can be 10 2 to 10 6 , and might partly com- 
pensate some of the decoherence effects due to the /i material 
itself. Increasing the SQUID dimension and decreasing the 
cavity volume will reduce the maximum allowed number of 
SQUIDs inside the cavity making it unadvantageous for quan- 
tum computing. However, one qubit is enough for the gen- 
eration of nonclassical photon states, our goal here. We note 
that Girvin etal. 1 17] proposed a different system in which the 
coupling of the photon-qubit can reach 10 8 Hz, corresponding 
to r c (0) - 10~ 9 s. We are considering how to generate non- 
classical photon states by using such a system. This scheme 
might not be easy to generalize in a straightforward manner to 
the flux qubit case. This because the interaction between the 
flux-qubit and the cavity field cannot be switched on and off in 
the same way for the charge qubit. However in some modified 
manner, it should be possible to generalize this scheme. 
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